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Abstract 

One of the main evolutionary stages of planet formation is the dynamical evolution of planetesimal 
disks. These disks are thought to evolve through gravitational encounters and physical collisions 
between single planetesimals. In recent years, many binary planetesimals have been observed in 
the Solar system, indicating that the binarity of planetesimals is high. However, current studies of 
planetesimal disks formation and evolution do not account for the role of binaries. Here we point out 
that gravitational encounters of binary planetesimals can have an important role in the evolution of 
planetesimal disks. Binary planetesimals catalyze close encounters between planetesimals, and can 
strongly enhance their collision rate. Binaries may also serve as additional heating source of the 
planetesimal disk, through the exchange of the binaries gravitational potential energy into the kinetic 
energy of planetesimals in the disk. 

Subject headings: planets and satellites: formation — protoplanetary disks — planets and satellites: 
dynamical evolution and stability — minor planets, asteroids: general 



1. INTRODUCTION 

Current models for the evolution of planetesimal disks 
consider three basic processes: viscous stirring, dynami- 
cal friction and coagulation or disruption through colli- 
sions (e.g. lSafronovll972l : [LissauerllT993l ; lGoldreich et al] 
[200l . These processes account for the rate at which 
single planetesimals encounter one another and collide, 
thereby changing the sizes (and masses) of planetesimals 
and their velocity distribution. None of these processes, 
however, takes into account the role of binary planetesi- 
mals. 

An important question when discussing the role 
of binaries is whether they exist at all in planetes- 
imal disks during their evolution. Evidently, in 
the Solar system many binary planet esimals (BPs) 
formed and survived to this day (see Gol dreich et al.l 
2001: 



WeidenschiUing 2002; Durda et al 



Funato et al... 2004: .Richardson fc Walsh,2006l:lLee et alJ 



2004 



i 



20071: ISchlichting fc Saril 120081 : IPerets fc Naozl 1200 
Nesvornv et al.l 120 lOL for possible formation and evo- 
lutionary scenarios). Large binary fractions (> 10 
%) are currently observed both for asteroids and 
trans-Neptu nian objects (TNO s) in the Solar s ystem 
(INoll et al.ll200a IKern fc El iotli 2006: M ann et al.l 120(171 : 

iWalsh. ,200^ 1. These are 



iRichardson fc WalshI |2006[ 



found to have separations of up to tens or hundreds 
times the radii of their single planetesimal components 
(for asteroids and TNOs, respectively; up to 0.1 the Hill 
radii for these objects). The currently observed large 
fractions of BPs in the Solar-system therefore serves as 
a basic motivation for our study. We assume that such 
BPs population exists in a given planetesimal disk, and 
study its role in the planetesimals evolution. We note 
that a variety of mechanisms we re suggested for the 
origin of BPs in the S olar sy stem ([Richardson fc WalshI 
120061 : iNesvornv et~all 120101 ). but only some limited 
aspects of their ev olution have been explored (e.g. 
IPerets fc Naod[2009l and references therein). 

|hperets@cfa.harvard.edu| 



The role of binaries in gravitational systems have been 
studie d in dep th in the context of stellar clusters ( |Hillsl 
119751 : iSpitzeil [19871 iHut et al.l 119921 ). In such systems 
stellar binaries are known to serve as important dynam- 
ical heating source of the clusters (so called binary heat- 
ing), which slows the collapse of the cluster. Binary stars 
are also known to have a major role in introducing stellar 
collisions and accelerate the growth of massive stars; it 
is thought that most stellar collisions are the result of 
encounters between binary and single (or other binary) 
stars. Here we suggest that BPs can play a similar role in 
planetesimal disks. The cross section for close encounters 
of binaries can be much larger than that of single plan- 
etesimals, thereby increasing their collision rates. The 
potential energy of close binaries could be exchanged into 
kinetic energy of the planetesimals (and vice- versa), and 
therefore serve as an additional heating (cooling) source 
for the planetesimal disk. In the following we explore 
both the process of binary induced collisions as well as 
binary heating. 

2. BINARY ENCOUNTERS 

The gravitational interactions between binary and sin- 
gle planetesimals can lead to various outcomes including 
binary disruptions, exchanges, resonant encounters and 
physical collisions. 

In binary-single encounters energy is exchanged be- 
tween the internal orbital energy of the binary, Etm-, and 
the kinetic energy of the incoming perturber, Ek- These 
are given by Ebin = Gmim2/'2a, where TOi,m2 are the 
binary-components masses and a is the binary mutual 
separation, and Ek — 0.5mw^, where m is the typical 
mass of the perturbing single planetesimals, and v is the 
relative velocity between the planetesimals (typically of 
the order of the velocity dispersion of the planetesimals). 
A binary is termed soft if \Ei,in\/Ek < 1 and hard if 
\E\/Ek > 1. On average, hard binaries get harder fol- 
lowing an encounter , and soft bin aries get softer (the so 
called Heggie's law: lHeggielll975f) . Soft binaries rapidly 
evaporate through such encounters and have a weak af- 
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feet on the energy budget of the system. However, en- 
counters with hard binaries lead, on average, to the sig- 
nificant loss of orbital energy from the binary, making 
the binary harder, following the scatter of the perturb- 
ing object into higher velocity. Therefore such a process 
leads to the increase of the planetesimals velocities and 
the dynamical heating of the system. Alternatively, a 
physical collision may occur between two or all of the 
planetesimals involved in the encounter. In the following 
we discuss the role of both these possibilities. 

2.1. Binary induced collisions 

The collision rate of a planetesimals is given by 

r — nav, (1) 

where n is the number density of planetesimals in the 
disk, a is the cross section for the collision and v is the 
relative velocity between the planetesimals. The cross 
section for a physical collision depends on the relative 
velocity between particles. When the relative velocity 
between planetesimals is much larger than the escape 
velocity of the most massive planetesimal participating 
in the encounter, gravitational focusing is negligible and 
the cross section is given by the projected cross sectional 
of the planetesimal a ttt^ (r is the planetesimal ra- 
dius). When the relative velocity is slower, gravitational 
focusing becomes dominant and we obtain the following 
cross section (cti) for two single planetesimal collisions 
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where v is the relative velocity and is the escape ve- 
locity from the planetesimal surface (with the right most 
term obtained for ^ u) . 

Binary-single encounters could involve complex tra- 
jectories, as the encounter now involves three bodies 
and the planetesimal trajectories can become chaotic 
(binary-binary encounters, not discussed here, interact 
in even more complex way). In such encounters the 
probability for a direct collision between any two (or 
even all) of the objects involved is high ly increased 
(jHills fc FullertonI [19801 : fFregeau et al.ll2004f ). and there- 
fore binaries serve as efficient catalysts for direct phys- 
ical collisions. The cross section for a physical colli- 
sion during a si ngle-binary en counters, CT2, is approxi- 
mately given bv iFregeau et al.l (2004; also Fregea.u, pri- 
vate communication a.nd (iSigurdsson fc PhinnevI 119931 : 
iValtonen fc Karttunenll2006h and references therein) 
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r/lO fcm 



a/2.14 X 103/cm 



0.65 
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where Vc is the critical velocity separating soft and hard 
binaries and a is the binary semi-major axis. Note that 
this equation was derived for stellar encounters, however, 
the dynamics of gravitational encounters is scale free, and 
can be scaled for the use of planetesimals mass objects. 
Although physical collisions also depend on the density of 
the objects, the typical average density of planetesimals 
(0.5 — 3 gr cm~^) is comparable to that of stars, and 
therefore scaling of mass and radii could be used. 

For simplicity we will consider equal mass planetesi- 
mals. For this case Vc is approximately the mutual or- 
bital velocity of the binary components (similar to the 



escape velocity from the binary at its separation, a; i.e. 
Vc (X We(f/a)^^^) and we therefore get 
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33 



0.35 



(4) 



As can be seen from this ratio, collisions during binary- 
single encounters could be tens up to hundreds of times 
more frequent than single-single encounters for typical 
BPs currently observed in the solar systerrQ. The ratio 
between collisions rates due to binary-single encounters 
(r2), vs. single-single encounters (Fi) is therefore given 

by 



r2 ^ n^a^v 

Fi TlifTiW 



39 



»2/»l 

0.3 



ajr 



0.35 



(5) 



where we took a binary fraction of — 712/ni — 0.3 in 
the normalization. 

Consequently, the rate of physical collisions and the 
mass growth of planetesimals (which is proportional to 

the collision rate, M oc F) could be boosted and domi- 
nated by binary-single encounters even for a low binary 
fraction. Moreover, binary-single encounters produce a 
different (power law) dependence of the growth rate on 
the planetesimals size. They also introduce dependencies 
on the binary parameters, namely the binaries separa- 
tion and binary fraction (which could be size dependent 
by themselves). The mass growth rate of planetesimals 
could therefore change both quantitatively (faster) and 
qualitatively due to existence of BPs. This could affect 
the size distribution of planetesimals in protoplanetary 
disks (or debris disks) and it s evolution. Henc e, the or- 
bital properties of binaries (jNaoz "eFaLllMot ). the size 
distribution of planetesimals, and their coupling could 
all be used, in principle, to characterize and constrain 
the evolution of the Solar system. 

Taking a conservative binary fraction of 10 percents for 
binary asteroids with typical separations of 10 times the 
planetesimal ra dius (e.g. typical of main belt asteroids; 
iRichardson fcW alsh 20Qi), we get r2/Fi - 7. The bi- 
nary fraction of TNO binaries could be higher than 30 
percents with typical separations of a ^ 50 r (e.g . typ- 
ical for binary TNOs; 'Rich ardson fc WalshI 120061 ). lead- 
ing to collision rate as high as ~ 39 times higher than 
the expected single-single physical collision rate. Given 
the likely higher binary fraction of planetesimals, both 
today and in the past, collision rates catalyzed by BPs 
during planet formation were likely to be even higher. 
Nevertheless, more detailed study of the formation and 
destruction mechanisms of BPs is required to asses this 
question quantitatively. 

We note that in the case of single-single collisions, the 
planetesimals are unbound prior to collision, whereas in 
the binary-single case the colliding planetesimals could 
have been marginally bound (e.g. during resonant en- 
counters) . The expected impact velocities during binary- 
single encounters are therefore likely to be lower, on aver- 
age, than those expected in single-single collisions. Lower 

^ Note that this ratio exceeds unity even as we approach to 
r = a; in these cases the flyby of a perturber can easily perturb 
the binary components into coUision, even if the perturber never 
crosses between the binary components. 
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velocity collisions are more likely to result in mass accre- 
tion rather than shattering of planetesimals (the likely 
outcome of higher velocity impacts), and are therefore 
more efficient for planetesimals growth. The accelerated 
collision rate may also imply that more collisions occur at 
earlier times, when the planetesimal disk is cooler, again 
leading to typically lower velocity impacts. 

2.2. Binary heating 

In the following we focus on the more energetically 
important hard-binaries. Encounters with hard binaries 
lead, on average, to the loss of orbital energy from the 
binary and the scattering of the perturbing object into 
higher velocity. There fore such a process lead s to the 
heati ng of the system (jHeggid ITOTa lHillj|1975t iSpitzeii 
Il987f ). In fact, in stellar clusters, this "binary heating" 
process is considered as one of the main processes gov- 
erning the evolution of the clusters (jHut et al.lll992[ ). We 
can now consider the effect of binary heating on the evo- 
lution of a planetesimal disk. We first note, however, 
that even the initial formation of binaries would heat 
the planetesimals disk (or the gas if the binaries formed 
due to gas interactions), since the binaries binding en- 
ergy have had to be transferred to the disk planetesimals 
upon their formation. 

2.2.1. Energy budget 

One can estimate the amount of potential energy reser- 
voir available for hard binaries. A binary can be hard- 
ened up to the point when it becomes a contact binary. 
The potential energy of such binary is of the order of 



Ebin 



2r 



(6) 



where r is the typical radius of a planetesimal of mass m. 
The energy extracted from a binary, which was initially 
at some typical separation a ^ r is of the order of 



^Eun — E 



fin - 



-E,, 



2r 



= E, 



bin • 



(7) 



If all this energy were to be gained by the planetesimals in 
the disk, they should have been excited to higher velocity 
dispersion, with 
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where Aw is the change in the velocity dispersion of the 
planetesimals, and N and N^m are the number of single 
and binary planetesimals, respectively. We therefore find 
that the change in velocity is 



At 



Nbinm _ Gm 



(9) 



where /f,i„ is the binary fraction. Comparing this to 
the Hill velocity, vh = [Gm/ RhY^^ (where Rh — 
A(3m/M0)^/'^ is the Hill radius at distance A from the 
Sun), we find that 



Av ~ \l fbin—VH- 

r 



(10) 



Typically, the radii of planetesimals are much smaller 
than the Hill radius by a few orders of magnitude. 
Therefore, even a small binary fraction potentially holds 
enough energy to heat up a disk to Hill and even super- 
Hill velocities. This could have important implications 
on the binary formation mechanism involved. For ex- 
ample, bi nary formation throug h dynamical friction dis- 
cussed bv lGoldreich et aTl ()2002|) , is efHcient only for rela- 
tively low velocity dispersions of the small planetesimals. 
This may imply that the high binary frequency observed 
in the Solar system is limited to large planetesimals; e.g. 
the energy extracted from the binaries formation is dis- 
tributed over a much larger mass of low mass planetesi- 
mals, which never formed binaries by themselves as these 
were all too soft to survive and contribute to the binary 
heating. Alternatively, binary formation could have been 
primordial, when they were embedded in gas, or during 
the collaps e of planetesimals swa rms forming single plan- 
etesimals (jNesvornv et al.l 120101 ). In these latter cases 
the extra energy from the binaries formation would be 
dissipated in the gas or radiated as heat during the grav- 
itational collapse. Better understanding of binaries for- 
mation is required for resolving this issue. 

2.2.2. Heating rate 

We now turn to the rate at which binaries heat the sys- 
tem. For simplicity we will consider a planetesimal disk 
composed of two types of objects: single planetesimals of 
mass TO, number density n, and typical velocity disper- 
sion V and hard binaries of mass rribin (where we assume 
the mass ratio, q, is typically high), number density ribin- 
The rate of energy gain b y the system from a siiigle hard 
binary is then given by (|Hills fc Fullertonl[T980l:TSpitieJ 
IT987l : lHm^lT99l 



dEh 



dt 



^D,,G^ 



''bin 
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where n is the number density of planetesimals, Venc is 
the typical encounter velocity (~ the velocity dispersion 
of the planetesimals), and D2 is some constant pre- factor, 
which will be discussed later on. We now use the relation 
dE/dt = Zmvdv / dt (the 3 factor comes from considering 
the velocity dispersion in a 3D system) to find expressions 
for the binary heating of the single planetesimals. Taking 



dE 
t— — 
dt 



dEb. 



-ribr. 



dt 
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and translating this into evolution of the velocity disper- 
sion (i.e. the energy extracted from the binaries is the 
kinetic energy gained by the planetesimals), we find 



rib, 

^dt'"'" 



1 dE, 



bin 



n 3mv dt 



^DoG 



2 ' '^binlTli^j^^ 
3vi^enc 



DoG 



2 nbinrribin 



(13) 

where numerical simulations (jHillsl [l"992| ) show that the 
constant pre-factor, D2, is of the order 6.7 for encounters 
of single mass objects (i.e. q 1 and nibin = 2m) and 
about twice as large for the case of more massive binaries 
(m mbin)- In the last term we assume the encounter 
velocity of of the same order of the velocity dispersion of 
the planetesimals. 
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We can compare the importance of binary heating with 
that of viscous heating by single planetesimals. Let us 
use a simple estimate for the viscous hea t ing, f ollow- 
ing a similar approach bv [Alexander et all ()2007l ). tak- 
ing a simplifying assumption that the velocity disper- 
sion is isotropic. This is not strictly valid, but it has 
been shown that the ratio of the radial and vertical 
velocity dispersions cannot become larger than 3 with- 
out the system becoming unsta ble ([Kulsrud et al.lllQTll : 
iPoliachenko fc Shukhmanlll977[ ) . The relaxation time for 
such a system is given by 



''relax 



CG^nm? In A 



(14) 



(e.g. iBinnev fc Tremaini Il987t iPapaloizou fc TerquemI 
[2OOII where n is the planetesimals density, v is the typical 
velocity dispersion of the planetesimals and In A is the 
Coulomb logarithm (where A ^ H/r\ H is the disk scale 
height) and C is an order-of-unity constant that depends 
on the ge ometry of the system, (fo r a spherical system 
C ~ 2.94. IBinnev fc Tremainill987[ ). Consequently, the 
relaxation of the system is governed by 
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with Di ~ 2C. Comparing with the binary heating term 
in Eq. [T2] we find both heating mechanism have similar 
dependencies on the disk and planetesimals properties. 
Taking the same velocity dispersion for both binaries and 
single planetesimals, and considering the same mass for 
all planetesimals (m;,m = 2to), the two heating terms 
differ only by some constant pre-factor, and the relative 
binary fraction, i.e. 



_ 4-P2 rtj,. 



V dt 



bin 
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Taken at face value AD2/D1 ^ 4.5, and therefore the bi- 
nary heating contribution could be comparable to that 
of viscous heating and dynamical friction, for realistic 
binary fractions of 0.1 — 0.2. Note that in the case of 
dynamical friction effect on low mass planetesimals, we 
should take the binary mass to be larger than the mass of 
the heated planetesimals and D2 would be twice as large 
(jHillsl [199^ . Given these uncertainties and the simpli- 
fications used in the above derivations, one should be 
cautious in using them at face value. These findings do, 
however, suggest that binary heating is likely not negligi- 
ble, but is also not likely to be the single most dominant 
mechanism for heating of the planetesimals disk. 

Binary heating would continue as long as hard bina- 
ries exist in the disk. However, as with the other forms 
of planetsimal disk heating, heating the disk puffs it up, 
resulting in a larger volume and hence lower number den- 
sity of planetesimals and lower encounter rates. 

In addition, as the velocity dispersion is increased, 
more binaries become soft, and could be disrupted. We 
note that all of the the suggested binary origin scenar- 
ios become less efficient at higher velocity dispersions, as 
the number density of planetesimals decrease and gravi- 
tational focusing becomes less effective. Currently, most 
of the observed BPs are soft in terms of encounters with 
similar size/mass planetesimals, and none of the binary 



planetesimals formation scenarios is currently effectiv e in 
producing new binaries (jRichardson fc Walshll2006t be- 
side the formation of low mass binaries through radiative 
spin up in the inner regions of the Solar system). There- 
fore binary formation and major collisional evolution had 
to proceed at early times, and the observed binary frac- 
tion today is representative only of the survived bina- 
ries, where as the earlier binary fraction was likely to be 
higher. 

3. CAVEATS 

Stellar vs. planetesimal encounters: The main 
caveat in our discussion of the role of BPs is the use 
of the stellar encounters approximation for estimating 
planetesimals encounters. Such approximations neglects 
the effect of the Sun on the encounters. 

The important regime for binary-single encounters is 
when the impact parameter of the encounter is smaller 
than the binary separation (and the Hill radius of the 
binary). Such encounters therefore occur mainly in 
the regime dominated by the mass of the planetesimals 
rather than the Sun. However, during resonant encoun- 
ters, planetesimals can be scattered to distances larger 
than the Hill radius, and still come back to re-encounter 
the binary in the absence of the Sun potential, at which 
point the tidal forces induced by the Sun could per- 
turb them. Most resonant orbits, however, are likely 
to be at smaller separations (see e.g. [Hills 1983, for a 
related, although different problem); in fact the prob- 
ability of being ejected to some large distance r ^ a 
is comparable to the prob ability of the binary becom- 
ing unbound altogether fsee lValtonen fc Karttunenll2006l 
chapter 8.3). Such perturbation and possible quenching 
of the interaction due to the Sun would therefore mostly 
affect wide binaries, with separations close to their Hill 
radius. Note, however, that when taking into account 
the gravitational pull of the Sun, the interaction even 
between single planetesimals could become resonant, ba- 
sically leading to the temporary captu re of the plan- 
etesimals and a long lived in teraction (jAstakhov et al.l 
[2001 iRafikov fc Slepianl[2OT0[ l. The additional interac- 
tion with the Sun can therefore increase the rate of reso- 
nant encounters in which the the BPs interact chaotically 
and therefore have a higher probability for physical col- 
lisions. Few body simulation of such binary-single plan- 
etesimals, which are beyond the scope of this letter could 
give a more quantitative picture of these processes. 

Collisional cooling: In the discussion of binary heat- 
ing we neglected the effect of physical collisions. The 
amount of energy damping through collisions is not very 
well understood for planetesimals and depend on the 
(unknown) coefficient of restitution of the planetesimals. 
The high collision rate during encounters could therefore 
serve an important role in the heating/cooling rate by 
binaries. We note that the importance of collisions is 
poorly known even in the much better studied literature 
of binaries in stellar clusters, where it was suggested to 
play an important role mostly for the hardest (tidally 
captured) binaries ( McMillan, 1986) making binary heat- 
ing much less efficient. 

4. DISCUSSION AND SUMMARY 

In this letter we pointed out the importance of BPs and 
their role in the evolution of protoplanetary and debris 
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disks. They could affect the evolution of such disks both 
through efficiently catalyzing physical collisions between 
planetesimals, and by serving as an additional planetes- 
inaal heating/cooling mechanism. 

Binary encounters provide a novel and efficient mecha- 
nism for planetesimals growth, suggesting a different de- 
pendence between the growth rate of planetesimals and 
their physical size. They also introduce dependencies on 
the binary properties, such as the binary fraction and 
semi-major axis (and their size dependent distribution). 
The evolution of the planetesimals size distribution in 
this case could therefore be qualitatively different than 
that envisioned in studies takin g into account only sin- 
gle p l anetesimals (e.K.[P ohnanvii l969t iDavis fc Farinellal 
Il997t iKenvon fc Bromlevll2004l and other related papers 
and references therein)]. 

The relative importance of binaries depends on the bi- 
nary fraction of planetesimals. Current observations of 
the Solar system show large BPs (asteroids; TNOs) popu- 
lations, even amongst the largest planetesimals/embryos, 
e.g. Pluto-Charon and the Earth-moon systems. BPs 
therefore exist even at late stages in the planetesimal 
disk, even up to the the scale of planetary embryos. BPs 
may therefore accelerate the mass growth of planetesi- 
mals both at early stages of evolution (when they were 



Alexander, R. D., Begelman, M. C, & Armitage, P. J. 2007, ApJ, 
654, 907 

Astakhov, S. A., Lee, E. A., & Farrelly, D. 2005, MNRAS, 360, 
401 

Binney, J. & Tremaine, S. 1987, Galactic dynamics 
Davis, D. R. & Farinella, P. 1997, Icarus, 125, 50 
Dohnanyi, J. S. 1969, JGR, 74, 2531 
Durda, D. D. et al. 2004, Icarus, 170, 243 

Fregeau, J. M., Cheung, P., Portegies Zwart, S. F., & Rasio, F. A. 

2004, MNRAS, 352, 1 
Funato, Y. et al. 2004, Nature, 427, 518 

Goldreich, P., Lithwick, Y., & Sari, R. 2002, Nature, 420, 643 

— . 2004, ARAA, 42, 549 

Heggie, D. C. 1975, MNRAS, 173, 729 

Hills, J. G. 1975, AJ, 80, 1075 

— . 1983, AJ, 88, 1857 

— . 1992, AJ, 103, 1955 

Hills, J. G. & FuUerton, L. W. 1980, AJ, 85, 1281 
Hut, P. et al. 1992, PASP, 104, 981 
Kenyon, S. J. & Bromley, B. C. 2004, AJ, 128, 1916 
Kern, S. D. & Elliot, J. L. 2006, Apjl, 643, L57 
Kulsrud, R. M., Mark, J. W. K., & Caruso, A. 1971, Ap&SS, 14, 
52 

Lee, E. A., Astakhov, S. A., & Farrelly, D. 2007, MNRAS, 379, 
229 



suggested to form; iGoldreich et al.ll2002[ ). and possibly 
even up to the stages relevant for the formation of the 
gas planets cores. 

We note that approaches used for the study of binaries 
in stellar clusters, could similarly be used to confront the 
new challenges raised by accounting for the role BPs, but 
we caution, and raise some caveats for their direct use in 
this context. 

Current studies of planet formation do not take into 
account BPs. Such additional component could be diffi- 
cult to account for in these studies (especially in simu- 
lations where following binary orbits is computationally 
expensive). Nevertheless, our findings suggest that not 
only that binaries are not negligible, but they may have 
an important role in the evolution of planetesimal disks, 
and therefore should not be ignored. We conclude that 
the inclusion of binaries in future studies could have im- 
portant implications for our understating of the evolution 
of planetesimal disks and planet formation. 
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